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ABSTRACT 

The  formation  of  hydraulic  jumps  in  a  rotating,    continuously- 
stratified,    Boussinesq  fluid  is  investigated.     A  scale  analysis  indicates 
that  jump-formation  cases  should  be  separated  from  non-jump  cases 
by  the  curve  F  -  ARo     where  Ro  and  F  are  the  initial  Rossby  and 
Froude  numbers,    respectively.      Numerical  solutions  verify  the  curve 
up  to  F  =  0.3  and  show  that  A  is  in  the  range  6.0  to  7.5.      The  solutions 
for  F  =  1.0  are  examined  in  considerable  detail.     The  results  of  this 
paper  extend  the  work  of  Houghton  and  Williams  and  Hori. 
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I.     INTRODUCTION 

The  formation  of  hydraulic  jumps  has  been  shown  for  a  one-layer, 
homogeneous  fluid  for  both  non- rotating  and  rotating  cases.      This 
experiment  will  consider  a  one -layer,    continuously  stratified  fluid  by 
first  investigating  jumps  in  a  non- rotating  coordinate  system  and  then 
extending  the  study  to  the  rotating  system. 

In  a  non- rotating  one -layer  fluid  a  hydraulic  jump  with  a  discon- 
tinuity in  the  velocity  can  arise  from  an  initial  state  having  any  scale. 
If  the  system  rotates,    there  may  be  a  limiting  initial  scale  beyond 
which  jumps   cannot  form  as   shown  by  Houghton  [1969]  and  Williams 
and  Hcri  [1970].      A.  similar   behavior  is  expected  ir:  a  continuous  ±v 
stratified  atmosphere.     In  non-dimensional  terms  this  corresponds 
to  a  critical  Rossby  number  (Ro  -  U/fJL,    where  U  is  the  initial  velocity 
scale,    L  is  the  initial  space  scale  and  f  is  the  Coriolis  parameter) 

below  which  jumps  are  excluded.     This  process  also  depends  on  the 

-1  h 

Froude  number  (F  =  U/[D(g0o  d0o/dz)    ],  where  D  is  the  vertical  scale 

and  0     is  the  mean  potential  temperature).     In  the  F-Ro  plane  there 

should  exist  a  critical  curve  that  would  separate  the  jump  cases  from 

the  non-jump  cases  as  has  been  shown  for  a  homogeneous  fluid  by 

Houghton  [1969]  and  Williams  and  Hori  [1970].     In  this  experiment  it 

was  shown  that  a  critical  curve  exists  and  is  valid  for  time  periods  of 

- 1  h 

order  L/[D(g0      dQQ/dz)    ].      For  longer  time  periods,    however,    jumps 


may  occur  in  the  non-jump  region,    indicating  that  the  critical  'curve' 
is  actually  a  zone  where  the  jump  time  increases  rapidly. 

There  are  several  possible  applications  of  this  work  in  meteorology 
and  oceanography.     Rotation  does  have  an  effect  on  mesoscale  motions 
in  the  atmosphere,   but  the  influence  is  not  dominant  as  it  is  in  large- 
scale  motions.     Tepper  [1952]  has  proposed  that  squall  lines  are 
modified  hydraulic  jumps.      Freeman  [1948],    Abdullah  [1966],    and 
Houghton  and  Kasahara  [1968]  have  proposed  other  mesoscale  hydraulic 
analogies.     Sasaki  [1959]  and  Ogura  and  Charney  [1962]  have  studied 
squall  lines  in  continuously  stratified  atmospheres. 

Internal  waves  in  the  ocean  are  affected  by  rotation,    and  the  break- 
ing cf  these  waves  should  be  described  by  this  model.     Rapidly  moving 
frontlike  temperature  disturbances  have  been  observed  in  the  ocean 
[Gowans,    1969]  and  these  may  be  internal  waves  modified  by  rotation. 

In  this  experiment  the  Boussinesq  equations  are  given  and  the 
modeling  relations  are  introduced  in  II.     In  the  model  the  time-dependent 
quantities  are  functions  of  x  and  z  only.      The  quantity  f  is  assumed  to 
be  constant  for  each  independent  case  study.      The  finite  difference 
approximations  are  presented  in  III.     A  scale  analysis  is  presented 
in  IV  which  is  useful  in  the  discussion  of  the  numerical  solutions.      The 
initial  conditions  are  given  in  V  and  the  numerical  solutions  are  pres- 
ented and  discussed  in  VI.      The  initial  state  consists  of  a  sinusoidal 
disturbance  (amplitude  dependent  on  F)  which  is  superimposed  on  the 
basic  state.     After  a  sufficient  period  of  time,    the  sine  wave  is  distorted 


in  such  a  way  that  a  hydraulic  jump  is  formed.     In  VII  conclusions  and 
recommendations  for  further  investigation  are  given. 


II.     BASIC  EQUATIONS 

In  this  experiment  the  Boussinesq  equations  are  employed  and 
the  domain  is  bounded  by  two  rigid  horizontal  planes.     The  compress- 
ibility of  the  atmosphere  is  neglected  in  the  Bossinesq  approximation. 
This  should  not  be  of  qualitative  importance  since  the  density  scale 
height  of  the  atmosphere  is  much  larger  than  the  vertical  scale  of  the 
typical  jump  zone.     The  hydrostatic  approximation  is  also  used  in 
this  experiment  and  it  should  be  good  as  long  as  the  horizontal  scale 
is  much  larger  than  the  vertical  scale  [Lamb,    1945,    p.    258]. 

The  Boussinesq  equations  which  were  obtained  by  Ogura  and 
T~':^11  }--i-  ^i^^^i  can  h*-  written  m  the' following  jum,  when  the  eartn s 


rotation  is  included: 

+  V'(vv)  +  17  (^V)  +  v<t>  +f  kxV»p,  (ii.  i) 


H 


|f  +  7-GV    +fe(w9)  =  0;  (n.2) 


\7-V  +lr  =0     and  (IL3) 


30  3_©  (II.  4) 


where: 

0O      =         constant  reference  potential  temperature, 

P        =         constant  reference  pressure, 
o 

k        =        R/Cp, 

6         =         T(PQ/P)      _  0O  =  departure  of  the  potential  temperature 
from  0O,    and 

0      =      cpe0(p/p0)*    +  gz. 

Heating  and  friction  have  been  neglected  in  these  equations  even  though 
they  may  be  important  under  certain  conditions. 

All  quantities  are  independent  of     y    in  this  model.      The  Coriolis 
parameter,    f,    is  assumed  to  be  constant  throughoiit  the  model  but  is 
varied  for  individual  cases.     A  Cartesian  coordinate  system  is  intro- 
duced with  x  increasing  toward  the  east  and  y  increasing  toward  the 
north.      The  x  and  y  velocity  components  are  u  and  v  respectively. 

For  this  model,    Equations  (II.  1 ) - (II.  3)  can  be  written  in  the 
following  form: 


at    +    i*  az 


iv       .     J{wv)    ,     ^(vw) 


+  £    -fv-O, 


«     +    *^v±  +  ^  +fw    =0  (IL6) 


H  ix  *Z 


ae       ijuG)      «w9).  0     and  (IL7) 

at   T    >x     T     « 

Ju      ^      iw  ^  (II.  8) 
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The  boundary  conditions  are 

w        =         0,    z  =  0,  D  (II.  9) 

where  D  is  the  distance  between  the  horizontal  plates. 

The  integral  of  Equation  (II.  4)  gives  the  departure  of  0  from  its 
vertical  average,  and  this  average  can  be  obtained  as  follows.  The 
following  averages  are  defined: 

where  2  n  L  is  the  periodicity  of  the  fields  in  x.      If  Equation  (II.  8)  is 
averaged  vertically  and  (II.  9)  employed 


><«> 


:0        •-     <«>=<(u)>.  (11.11) 


Similarly,    the  vertical  average  of  (II.  5)  is 

K")    .    K"u>     .     ><0>        rn  „  Q  (II.  12) 

it  *x  ^X  ^V>    " 

If  this  equation  is  in  turn  averaged  with  x  and  if  periodicity  in  x  is 
imposed,    Equation  (II.  12)  can  be  written 

^    -  f <tW>  =  0.  in- 13) 

The  quantity  l(V)l  represents  the  net  mass  flux  in  the  y  direction.     If 
/(V)^    were  different  from  zero,  /(y— J)  could  not  be  zero  if  a  wall 
were  present  at  some  value  of  y.     Since  this  would  violate  the  boundary 
conditions  [Equation  (II.  9)],    the  following  is  a  requirement  in  this 
experiment: 
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<(v)>   r  0. 


(II.  14) 


It  then  follows  directly  from  (II.  13)  that 


d<(u)> 


=  0. 


(11.15) 


Now  from  (II.  11) 
reduces  (II.  12)  to 


><u>         cj<(u)> 


d* 


which  together  with  (II.  15) 
+    £<v>.  (II.  16) 


The  corresponding  equation  is  obtained  by  averaging  Equation  (II.  6)  in 
x  and  z  and  by  introducing  (II.  14),    i.  e.  , 

*<(*>>*  o.  (II-17) 

Equation  (II.  17)  shows  that  the  x-  and  z-averaged  zonal  wind  is  zero. 
if  Equations  (II.  4),    (II.  16)  and  (II    L7)  are  used,    tnua-or^  (II.  5)  and 
(II.  6)  can  be  rewritten  in  the  forms 


it         >x  ax  ;i  7Vk      x  /y 


+  ±-  (tf  -<*>)  =  o    and 


a* 


(II.  18) 


(II.  19) 


By  integrating  the  hydrostatic  equation  [Equation  (II.  4)]  the  following 
expression  for    $"(9)  is  obtained: 


*-<#>«  i[fedz  -(f}^)\ 


(11.20) 


The  following  energy  integral  is  an  invariant  for  this  model: 


e-  [<(«*)>  *<^)]A    -  J  <(?))• 


(11.21) 
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III.      FINITE  DIFFERENCE  APPROXIMATIONS 

Equations   (II.  7),    (II.  8),    (II.  18)  and  (II.  19)  were  solved  numerically 
by  introducing  finite  differences  in  x,    z  and  t.     The  solutions  were 
obtained  by  a  simple  marching  process,   but  if  the  horizontal  variations 
were  two-dimensional,    Poisson's  equation  would  have  to  be  solved  at 
each  time  step  [Smagorinski,    1958;  Hinkelmann,    1959]. 

The  vertical  domain  was  divided  into  2K  layers  of  equal  thickness, 
and  the  surfaces  separating  the  layers  were  denoted  by  k.      The  arrange, 
ment  of  variables  on  the  bounding  surfaces  is  shown  in  Figure   1.     The 
divergence  of  the  vertical  flux  of  the  scalar  S  at  z  =  z,    is  approximated 


y  s 


XwS) 


(III.  1) 


LK    ////////////////////  T,D 


k-.K-l 


w 


W 


k  =  |    w 


k»4  "v.0^ 


llllllllll 1 1 1  III  III 


Fig.    1.     Arrangement  of  variables  in  z. 
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The  finite  difference  approximation  to  the  integral  in  Equation  (11.20) 
used  is 

>-i   L 


f\dz=I:''l^M(z.ti.Zj^cj  ^i.        (IIL2) 

The  value  of  C  =  J  *Q4l  is  immaterial  since  it  will  disappear  when  the 


vertical  average  is  removed.      The  finite  difference  integral  of  Equation 
(II.  8)  used  is 

The  vertical  grid  structure  was  introduced  by  Lorenz  [i960]  for  use  in 
the  balance  equations. 

All  of  the  dependent  variables  occur  at  each  x  grid  point.      The  x 
derivatives  were  approximated  as 

(**'*  2Ax  am*  (IIL4) 


[ 


(U,>|-*"U/)(S,>lfS;)  -(u,-  +  U,.{)fS,fS,-l) 

4.  Ax  •     (IIL5) 


Here  i  denotes  the  grid  index  and  Ax  the  mesh  length.      This  is  a  variant 
of  the  finite  difference  scheme  that  Arakawa  [1966]  developed  [see  also 
Lilly,    1965].     The  Arakawa  scheme  is  designed  to  avoid  the  non-linear 
computational  instability  which  was  discovered  by  Phillips  [1959].     In 
particular,    the  finite  difference  approximation  to  the  integral, 
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is  zero  when  (III.  1)  and  (III.  5)  are  used.      The  integral  I  is  always 
zero  in  the  continuous  case.     If  the  domain  averages  of  u    ,    v     and  9^ 
are  conserved  by  the  non-linear  terms  throughout  a  numerical  integra- 
tion,   then  no  non- linear  computational  instability  can  occur.     When 
finite  differences  in  time  are  employed,    these  quantities  will  not  be 
identically  conserved  and  instability  may  arise. 

The  continuous  equations  conserve  the  total  energy  E,    which  is 
given  by  Equation  (11.21)  ,     If  non-linear  computational  instability 
occurs,    E  cannot  be  conserved.     It  was  found  that  E  remained  nearly 
constant  in  the  numerical  solutions  in  this  experiment,    which  indicates 
that  the  spurious  computational  production  of  energy  was  small.      This 
was  important  because  it  was  desired  that  the  computational  scheme 
would  not  introduce  a  small-scale  energy  source  which  could  produce 
a  false  jump. 

Centered  time  differences  were  used  throughout  this  experiment 
with  the  exception  of  the  first  time  step  which  was  a  forward  time  step. 
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IV.     SCALE  ANALYSIS 


A  scale  analysis  is  now  performed  that  will  indicate  the  dependence 
of  the  solutions  on  the  Rossby  number  and  Froude  number.      The  inde- 
pendent variables  are  written 

t  =  Tt',    x  =  Lx'  and  z  =  Dz\  (IV.  1) 

The  dependent  variables  are  scaled  as  follows: 

u  =  Uu',    v  =  Vv',     w=Ww\       0  =  UDN0  '  and 

6  =  (UN9o/g)0>  .  (IV.  2) 

The  Rossby  number  and  Froude  number  are  defined  as 

Ro  =  U/fL    and  F    =     U/DN  (IV.  3) 


ivr  ;„    j  ..  r: 


N  = 


"z 


[lie] 


(IV.  4) 


The  time  scale  T  and  the  v-scale  V  are  functions  of  Ro  and  F  while 
the  w-scale  is  W  =  UD/L. 

An  examination  of  the  case  where  F£  1  and  F&  Ro  is  now  given. 
This   should  parallel  the  results  obtained  by  Houghton  [1969].      The 
appropriate  time  scale  is  T  =  L/DN,    and  the  v_scale  is  V  =  (F/Ro)U. 
The  later  was  obtained  following  Houghton  by  balancing     yr       with  fu 
in  Equation  (II.  6).      The  non-dimensional  equations  for  this  case  are 


*«'  r- 


>("'«')  i(u'w') 


lY' 


iz' 


*Y' 


r2 


(IV.  5) 
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*v' 

77,   +  F 

at 


i(u'v')  ^(v'v//) 


>*' 


*Z' 


+     U    *   0 


(IV.  6) 


if 


+  F 


+ 


>*' 


az' 


t  w'  =  0, 


(IV.  7) 


-  +  jp  -o     wd 


*z' 


(IV.  8) 


(IV.  9) 


Hydraulic  jumps  are  formed  through  the  action  of  the  non-linear  terms 

.    ___   1--7   c\     .'TT..r    '.  \  _  .  ..-i  ,' i  v    '.'  \       hsjcv.  when  these  terms  ca.j.c 
small  (F«    1),    they  can  produce  a  jump  in  a  non- rotating  system.     If 
a  hydraulic  jump  is  to  be  prevented,    the  Coriolis   term  in  (IV.  5)  must 
dominate  the  non-linear  advective  terms,   which  implies  that 


(IV.  10) 
The  curve  dividing  the  jump  region  from  the  no -jump  region  should 


r2 


then  have  the  form 


F  =  A  R.a 


(IV.  11) 


where  the  constant  A  should  be  appreciably  larger  than  1.     Houghton's 
analytical  curve  corresponds  to  A  ^    6.5;  Williams  and  Hori  [1970] 
found  that  A  =  6.  0  to  7.  5  was  in  agreement  with  their  numerical 
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solutions.  As  in  Williams1  and  Hori's  experiment,  with  A  in  the  range 
6.0  to  7.5,  the  above  scale  analysis  does  not  apply  for  F  r\j  1,  since  F 
is  then  greater  than  Ro. 

Two  regions  in  Ro  are  of  interest  when  F^v    1.     If  Ro  rv   1,   all 
terms  in  Equations  (IV.  5),    (IV.  6)  and  (IV.  7)  are  of  the  same  order 
and  a  jump  would  be  expected  to  form.     However,   when  Ro((    1,    the 
appropriate  time  scale  is  T  =  f-1  and  the  v- scale  is  V  =  U.      The  non- 
dimensional  equations  then  become 


>u'  n 


ax'  dz'  f    dx' 


-  V    *0, 


+  u'  r  o     and 


^tR.[— +  -TF-  +7WJ    =° 


(IV.  12) 


(IV.  13) 


(IV.  14) 


along  with  Equations  (IV.  8)  and  (IV.  9).     If  the  terms  in  Ro  are  neglected, 
the  resulting  equations  describe  an  inertial  oscillation  in  u  and  v;  the 
temperature  field  is  independent  of  time. 

The  scale  analysis  indicates  that,   when  F  is  small,    the  jump 

2 
boundary  is  given  by  F  =  ARo    ,    and  that  as   F  approaches   1,    the 

boundary  will  be  between  Ro  =1.0  and  Ro  =  0.  1.      These  resxilts  were 

tested  by  direct  numerical  integration  of  the  basic  equations. 


V.     INITIAL  CONDITIONS 

The  initial  x-averaged  stability   r—    is  independent  of  z.     A  wave 
sinusoidal  in  x,    the  amplitude  of  which  is  directly  dependent  on  the 
Froude  number,    is  superimposed  upon  this  basic  state. 

The  initial  fields  for  this  experiment  are  given  by  the  equations 

u(«,0)  =  U  cosiki)  {cos[£  (z  -§)]  - 1} , 

v(«,0)=0    and  (V'I) 

e^W.-^Sr  e«(l.»)  $in[|(2-fj-g(z-f). 

The  constants  which  define  the  initial  state  have  the  following 
numerical  values: 

0  *  3"  J<w,,; 

ff  =  4.0  °K  lew!1    and  (v.z) 

1  -.  0.03269   m  scc*20^1. 

The  wavelength  of  the  disturbance,    L,    is  equal  to   1,  000  km. 

After  finite  differences  are  introduced,    it  is  necessary  to  adjust 
the  values  of  u  which  are  given  by  (V.  1)  in  order  that 

This  merely  states  that  the  finite  difference  vertical  velocity  at  the 


upper  boundary  must  be  zero. 
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The  following  initial  conditions  were  also  examined: 

uh.oNUcojMcosPr). 

v(x,o)  s  o    and  <v-4> 

eM)..?!pie«(fc,)$i»(f)  +  g(f.j). 

When  f  =  0  they  lead  to  an  exact  finite  amplitude  wave  solution 
which  moves  with  the  linear  gravity  wave  speed.     In  this  solution  the 
non-linear  terms  in  Equations  (II.  7)  and  (II.  18)  cancel  out  identically. 
In  this  case  no  jumps  can  occur  although  they  were  observed  for  f  £  0, 
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VI.     NUMERICAL  SOLUTIONS 

The  numerical  solutions  are  presented  in  non-dimensional  form. 
The  time  for  all  results  presented  is  given  in  units  of  the  time  scale, 
L/DN.     Most  of  the  results  are  for  F/v    1.      Enough  results  are  pre- 
sented for  F  <(    1  in  order  to  determine  the  jump  curve.     This  will 
provide  a  basis  for  comparison  with  the  work  done  previously  by 
Houghton  [1969]  and  Williams  and  Hori  [1970]  in  the  one-layer,    homo- 
geneous system. 

A  quantity  xiseful  in  the  analysis  of  the  jump  criterion  is  defined 


as 


The  quantity  r-r    is  approximated  by  (u^+j  -  u^)/  a  x.     Here  the  sub- 
scripts max  and  min  refer  to  the  maximum  and  minimum  overall 
points  at  a  given  time.      The  non-dimensional  quantity  d  is  a  measure 
of  the  width  of  the  zone  through  which  u  varies  most  rapidly.      It  will 
be  an  overestimate  of  this  width  when  the  values  of  u  at  the  edges  of 

the  zone  are  not  near  u  and  u     .  However,    if  a  discontinuity 

max  mm.  '  ' 

such  as  a  pressure  jump  does  occur,    d  must  approach  zero,    since  it 

is  inversely  proportional  to  the  gradient  of  u. 

Figure  2  shows  the  variation  of  d  as  a  function  of  time  (time  unit 

L/DN)  for  F  =  1.  0  and  Ro  =  00.     Results  are  shown  for  three  values  of 

A  x  =  2  7v  L/n  corresponding  to  n  -  20,    40  and  80.     As  can  be  seen,    d 

decreases  to  a  minimum  at  approximately  DNt/L  =  0.4. 
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Detailed  solutions  are  presented  for  three  cases  in  which  F  =  1. 
In  these  numerical  solutions       Ax  -  2-  L/n,   where  n  =  80.      Figure  3 
shows  the  time  evolution  of  the  non-dimensional  variables  u'  -  u/U 
v'  =  v/U  and  0'  =  (0  _  ~Q)  /(UN0    /g)  [where  "0  represents  the  horizontal 
average]  for  Ro  =  1.     The  u  and  0  fields  show  that  a  hydraulic  jump  is 
forming  at  DNt/L  =  0.  37.     The  v  field,    v/hich  begins  at  zero,    grows  in 
opposite  phase  to  the  u  and  0  fields  through  the  action  of  the  Coriolis 
force.     Once  the  v  field  is  present,    the  Coriolis  force  will  tend  to 
reduce  the  amplitude  of  the  u  field  while  the  amplitude  of  the  0  field 
remains  relatively  unchanged  throughout  the  period  of  integration.      For 
this  value  of  the    Rossby  number  the  Coriolis  force  was  unable  to 
prevent  the  formation  of  the  hydraulic  jump. 

Figure  4  contains  the  evolution  of  the  non-dimensional  variables 
u1,   v'  and  0'  for  Ro  =  0.  1.     It  appears  that  an  inertial  oscillation  is 
present  with  the  0  field  remaining  very  nearly  constant  throughout  the 
integration  period.     It  was  found  that  when  the  integration  time  was 
extended,    though,    the  fields  never  were  restored  to  the  original  as 
was  shown  for  the  one-layer,    homogeneous  case  by  Williams  and 
Hori  [1970]. 

The  results  for  Ro  =  0.  35  are  presented  in  Figure  5  using  the 
same  format  as  above.     In  this  case  there  is  no  evidence  of  any 
reduction  in  the  u  field  as  the  jump  time  is  approached. 
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Fig.    3.      The  fields  u',    v1  and 

9'  at  various  times  for  F=1.0 

and  Ro  =  l.  0. 
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Fig.    4.      The  fields  u\    v'  and 

9'  at  various  times  for   F=l.  0 

and  Ro  =  0.  1. 


2 

Fig.    5.      The  fields  u',    v'  and 
9'  at  various  times  for   F=1.0 
and  Ro=0. 35. 
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Results  for  F  =  0.  04  and  R.o  =  0.  09  are  contained  in  Figure  6 
showing  the  evolution  of  u',    v'  and  9'  for  a  typical  jump  case  when 
F  «  1.      Figure  7  contains  the  evolution  of  the  variables  for  F  =  0.  04 
and  Ro  =  0.  07  giving  the  results  for  a  non-jump  case  when  F  is  again 
much  less  than  1. 

In  Figure  8  and  Figure  9  the  results  are  presented  for  F  -  0.  3 
and  Ro  =  0.3  and  0.  1,    respectively,    showing  the  time  evolution  of 
u',    v'  and  0'  when  F  is  in  the  range  0.  1  to  1.  0. 

The  above  figures  are  taken  from  the  results  of  the  lowest  layer 
in  the  model.      Figure  10  is  included  in  order  to  show  what  is  taking 
place  in  the  other  layers  at  time  zero  and  at  an  approximate  jump 
time.      The  case  presented  is  for   F  =   1.  0  and  Ro  =   1.  0.     As  can  be 
seen,    only  the  lower  half  of  the  model  is  shown.      This  provides  a 
picture  of  what  is  taking  place  in  the  most  active  of  the  layers  in  the 
model. 

The  d  procedure  was  used  to  determine  the  curve  separating  the 
jump  cases  from  the  non-jump  cases  over  the  range  F  =  0.  04  up  to 
F  =  1.0.      These  results  are  shown  in  Figure   11,    with  the  logarithm 

of  F  plotted  against  the  logarithm  of  Ro.      The  lines   F  =  6.  ORo     and 

2 

F  =  7.  5Ro     are  included  for  convenience.      These  results  for  Equation 

(IV.  11),   which  was  derived  by  scale  analysis  leading  to  the  conclusion 
that  A  is  in  the  range  6.  0  to  7.  5.      This  relation  does  not  hold  for  F^v    1, 
but  this  was  to  be  expected  since  the  analysis  reqtiired  that  F  £    Ro. 
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Fig.    6.      The  fields  u',    v'  and 
0'  at  various   ames  for  F=0.  04 
and  R-j^O.  09. 
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Fig.    7.      The  fields  u',    v'  and 
9  at  various  times  for  F=0. 04 
and  Ro  =  0.  07. 
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Fig.    8.      The  fields  u',    v1  and 
0'  at  various  times  for   F=0.  3 
and  Ro=0.  3. 
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Fig.    9.      The  fields  u',    v'  and 
G'  at  various  times  for   F=0.  3 
and  Ro  =  0.  1. 
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Fig.    10.      The  lowest  five  layers  of  the  model  at  time  zero  and  at 
approximate  jump  time  for   F=l.  0  and  Ro=l.  0. 
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Fig.    11.     Summary  of  jump-formation  experiments  with  the  logarithm  of 
F  plotted  against  the  logarithm  of  Ro.      The  crosses   indicate  jump 
occurrence  and  the  dots  indicate  cases  where  no  jumps  were  found. 
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The  non-jump  solutions  included  in  Figure  11  are  determined  over 
an  integration  time  of  L/DN. 
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VII.      CONCLUSIONS  AND  RECOMMENDATIONS 

The  formation  of  hydraulic  jumps  in  a  rotating,    continuously- 
stratified  fluid  has  been  investigated.     A  scale  analysis  indicated  that 
jump  cases  should  be  separated  from  non-jump  cases  by  the  curve 
F  =  ARo     in  the  F-Ro  plane.     The  Rossby  number  Ro  and  the  Froude 
number  F  were  utilized  as  direct  inputs  to  the  equations  and  used  as 
a  means  of  determining  the  initial  scales  (U  =  FDN).     It  was  determined 
that  the  results  of  this  experiment  (with  F  appropriately  redefined) 
closely  paralleled  the  results  of  Houghton  [1969]  and  Williams  and 
Hori  [1970]  in  which  one-layer,    homogeneous  fluids  were  studied. 
The  numerical  solutions    'jf  this  experiment  imply  that  -m  is  in  the 
range  6.  0-7.  5  and  that  the  relation  F  =  ARo     holds  up  to  at  least 
F  =  0.3.     Extrapolation  to  F  =  1.0  showed  jumps  occurring  for 
smaller  values  of  Ro  than  were  given  by  the  above  curve.      Thus  it  can 
be  seen  that  the  jump  criterion  becomes  more  independent  of  Fas 
F  =  1.  0  is  approached. 

This  work  should  be  extended  to  an  atmosphere  in  which  the  jump 
region  is  bounded  by  a  stable  layer  rather  than  a  rigid  plane  as  in 
this  case.     Such  a  study  could  require  additional  layers  with  further 
consideration  given  to  layers  of  varying  thickness.      Further  considera- 
tion should  be  given  to  moisture  in  future  experiments. 
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